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We study incompressible ground states of bosons in a two-dimensional rotating square optical 
lattice. The system can be described by the Bose-Hubbard model in an effective uniform magnetic 
field present due to the lattice rotation. To study ground states of the system, we map it to a 
frustrated spin model, followed by Schwinger boson mean field theory and projective symmetry group 
analysis. Using symmetry analysis we identify bosonic fractional quantum Hall states, predicted for 
bosonic atoms in rotating optical lattices, with possible stable gapped spin liquid states within the 
Schwinger boson formalism. In particular, we find that previously found fractional quantum Hall 
states induced by the lattice potential, and with no counterpart in the continuum [G. MoUer, and 
N. R. Cooper, Phys. Rev. Lett. 103, 105303 (2009)], correspond to "vr flux" spin liquid states of 
the frustrated spin model. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices have provided new 
opportunities to experimentally realise and study a wide 
range of complex quantum models under controllable 
conditions ji"— Many studies of atomic systems in artificial 
gauge potentials were motivated by the possibility to re- 
alize strongly correlated fractional quantum Hall (FQH) 
states4"J^ On the other hand, quantum Hall effects^^"^* 
have attracted much attention since their discovery, and 
are still a very active field of study. In addition to being 
intrinsically interesting systems, they are also being in- 
vestigated from the point of view of topological quantum 
computing: excitations in some FQH systems exhibit 
fractional, anyonic statistics and can potentially perform 
fault-tolerant quantum computationJ^iSS 

An artificial gauge field for neutral atoms can be cre- 
ated, for example, by rotating the latticei^ i^'^'^^" — The 
rotation of the system is equivalent, in the corotating 
frame, to the introduction of an effective magnetic field 
proportional to the rotation frequencyi^i^ However, the 
system enters the FQH regime at a very high rotation 
frequency or a low atomic density which are hard to 
achieve experimentally. This has motivated development 
of many alternative scheme o^'^'^"^" — to introduce a vec- 
tor potential to a system of neutral atoms, for example, 
using Raman assisted hopping combined with lattice ac- 
celeration or an inhomogeneous static electric field, or 
using oscillating quadrupole fields.— 

Prior to studies of FQH states in optical lattices, such 
states have been predicted for rotating Bose-Einstein 
condensates in harmonic magnetic traps The rota- 
tion (effective magnetic field) introduces vortices into 
the condensate and eventually, at sufficiently high rota- 
tion frequency, leads to formation of an Abrikosov vortex 
lattice.'*''^'' When the number of vortices becomes compa- 
rable to the number of bosons the system can enter into 
a FQH state. However, since the interactions among the 
bosonic atoms in the magnetic traps are weak, the FQH 



ground state is separated from the excited states by a 
small energy gap. On the other hand, in optical lattices 
interaction energies among the atoms are much larger due 
to confinement of atoms in a smaller volume. This leads 
to more robust FQH states separated from the excited 
states by much larger energy gaps than for the atoms in 
harmonic traps. In addition, the presence of the lattice 
leads to new interesting physical effects. In particular, 
new FQH states with no counterpart in the continuum 
limitiS, were predicted for systems of bosons in the pres- 
ence of the lattice potential and effective magnetic field. 

In the presence of a tight-binding lattice and a uniform 
magnetic field the single-particle energy levels change 
from simple Landau levels in the absence of the lat- 
tice into the Hofstadter butterfly spectrumi^i^ The en- 
ergy spectrum depends on the magnetic flux a per lat- 
tice plaquette (measured in units of the flux quantum 
4>o = 27rh/Q, where Q is the effective charge of the par- 
ticle). 

For rational a — p/q (with p and q coprime) the spec- 
trum splits into q bands and each state is q fold degener- 
ate. When the single-particle energies are plotted against 
the flux a the resulting Hofstadter butterfly has a fractal 
structure. For a ^ 1 the low-lying bands become the 
Landau levels of the continuum. For bosonic atoms with 
repulsive interactions, in the presence of a lattice and a 
uniform effective magnetic field with a <C 1 one thus ex- 
pects to find the same states as in a continuum bosonic 
quantum Hall system. Some of the states that appear in 
the continuum are, for example, bosonic v = 1/2 Laugh- 
lin state and ly = 1, 3/2, ... Read-Rezayi states^">2§ with 
1^ = n/a and n being the average particle density per 
lattice site. It has been shown that at sufficiently low 
particle density the lattice has a negligible effect on the 
nature of the continuum ly — 1/2 statei^iii However, in 
the presence of the lattice additional novel FQH states 
with no counterpart in the continuum limit have also 
been predicted 

In the present paper we study incompressible ground 
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states of interacting bosons in a two dimensional rotating 
square optical lattice. Bosonic atoms in a deep optical 
lattice can be described by a Bose-Hubbard model. An 
effective vector potential, due to the lattice rotation, in- 
troduces an Aharonov-Bohm phase for the bosons hop- 
ping from site to site and the complex tunelling cou- 
plings appear in the Hubbard Hamiltonian. In the regime 
of strong on-site interaction the Bose-Hubbard Hamilto- 
nian can be mapped to an effective frustrated spin model. 
To study the spin model, we use Schwinger boson mean 
field theory^"— and projective symmetry group (PSG) 
analysis4i"— 

Schwinger boson mean-field theory allows description 
of both ordered and disordered phases, that is, the result- 
ing mean field Hamiltonian does not have any preferred 
direction in spin space. Magnetic ordering is identified 
as Bose-condensation of Schwinger bosons. For such con- 
densed phases the energy excitation spectrum is gapless, 
while a gapped spectrum corresponds to incompressible 
spin liquid phases. In this work, we focus on the possible 
noncondensed ground states of the system, which do not 
break any symmetries. For such states, PSG anlysis can 
be used to determine different mean field ansatze that 
reflect all the physical symmetries of the system and to 
classify possible spin liquid states. 

Unlike conventional states which can be distinguished 
by patterns of broken symmetry, different spin liquid 
phases that are completely symmetric can be distin- 
guished by considering how the symmetries are realized 
in those phases. Such symmetric spin liquid phases can 
be classified using PSG analysis. The main idea of the 
PSG analysis of the mean field states in the Schwinger bo- 
son theory is that the mean field ansatz preserves all the 
symmetries of the spin Hamiltonian when it is invariant 
under transformations that are combined physical sym- 
metry and U{1) gauge group transformations. In other 
words, symmetries of the spin model can be preserved in 
the mean field state in some indirect way and invariance 
of the mean field ansatz under all symmetry transforma- 
tions is a sufficient, however not necessary condition. The 
set of all transformations that leave a mean field ansatz 
invariant is called the PSG. 

Different allowed sets of combined transformations can 
be found by considering algebraic relations among sym- 
metry group elements. Those different allowed sets of 
combined transformations then characterize physically 
distinct symmetric spin states. The bosonic FQH states 
of the system can then be identified with different spin 
liquid states of the effective spin model. In particular, 
our results show that lattice induced FQH states of the 
system correspond to "tt flux" spin liquid states of the 
effective spin model. 

The paper is organized as follows. In Sec. |ll]we intro- 
duce the Bose-Hubbard model for the system and outline 
the mapping to the quantum spin model. In Sec. IIIII we 
describe Schwinger boson mean fleld theory for the effec- 
tive quantum spin model. The magnetic symmetry group 
operators that correspond to the symmetries of the sys- 



tem are defined in Sec. IIVI In Sec. |V] we find possible 
noncondensed spin liquid ground states of the system us- 
ing PSG analysis. Relation between the spin liquid states 
found in Sec. IVland bosonic FQH states of the system is 
investigated in Sec. IVIl In the final section, we draw our 
conclusions. 



II. EFFECTIVE FRUSTRATED SPIN 
HAMILTONIAN 

We consider the system of bosonic atoms in a two- 
dimensional optical lattice and in the presence of an ef- 
fective vector potential introduced, for example, through 
lattice rotation. In the limit of weak tunneling t between 
wells within the lattice, compared to the level spacing in 
each well, the system can be described by a single-band 
Bose-Hubbard model on a square lattice with a complex 
hopping matrix element: 

i^Hubbard - + V (1) 

with 

i i 

V ^ -tJ2 {e"^''a]a, + e-'^-a|aj) , (2) 

where Oi and al are bosonic field operators on site i, U 
is the repulsive energy of two atoms in one well, /i is the 
chemical potential and (ij) denotes nearest-neighbor sites 
i and j. The complex tunneling couplings appear in the 
Hubbard Hamiltonian due to the presence of the effective 
vector potential A that introduces an Aharonov-Bohm 
phase for the boson hopping from site to site. When an 
atom moves from a lattice site at ri to a neighbouring 
site at rj , it will gain an Aharonov-Bohm phase: 

= r A-df. (3) 

For a rotating lattice A — mil x r/h, where (l is the ro- 
tation frequency and m is the mass of the atom. We con- 
sider the case of the uniform magnetic field B = V x A = 
Bz. The vector potential introduces frustration in the 
atomic motion if the phase twists around each plaquette 
add p to 27ra for some noninteger a. The frustration pa- 
rameter a is defined as the flux per plaquette in units of 
27r: 

a= _ / B-d5plaq= — ^</)^J, (4) 

zTT J zvr ^ — ' 

plaq 

where the integration is over the surface of a lattice pla- 
quette and the sum is performed anticlockwise over the 
edges of the square plaquette. Due to the periodicity un- 
der a — a-f 1 we can restrict the values of a to < a < 1. 
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The frustration is maximal at a = 1/2. We choose the 
Landau gauge A = B {0, x, 0) so that the Aharonov- 
Bohm phase (pij is zero on all horizontal bonds of the 
lattice, (l){x,y){x±i,v) = 0, and (l){x,y)(x,y±i) = ±2TTax. 
We further focus on the limit of weak tunneling com- 
pared to the repulsive energy U for two atoms in one well, 
t <^ U , where strong interparticle repulsion can lead to 
strongly correlated ground states. In such a limit the site 
occupation can be restricted to zero and one boson and 
the Bose-Hubbard Hamiltonian (jlj can be mapped onto 
a spin-1/2 frustrated XY model. The two 5*2 states of 
the pseudospin at a lattice site i correspond to whether a 
lattice contains a boson or not. The strongly correlated 
phases of the system that we study in the present paper 
will correspond to spin-liquid phases of the effective frus- 
trated quantum spin model. 

The mapping of the Hubbard Hamiltonian ([T]) in the limit 
oit -^U to an effective spin Hamiltonian is possible be- 
cause hard-core boson operators have the same commu- 
tation relations as spin-1/2 operators. The operators on 
different sites commute and operators on the same site 
anticommute. The spin raising and lowering operators 
correspond to the creation and annihilation operators of 
hard-core bosons and the motion of the atoms translates 
to pseudospin exchange. The effective spin Hamiltonian 
is: 

(ij) i 

(5) 

where J = 2t, = Sf±iSf are spin-1/2 operators, and 
h — ^ represents an effective Zeeman field. Note that this 
is a ferromagnet in the absence of frustration (0^^ = 0) . 
Also, it can be easily verified that the Hamiltonian ^ 
has local gauge invariance. If we change the gauge, 
A —i' A + Vx, then the Hamiltonian stays unchanged 
if the boson and spin operators pick up a phase change, 
^ e'(x^-»V^J , a^ e'^'a^ and S'r -> e'»5r. The 
change of gauge corresponds to a spin rotation of Xi in 
the xy plane. 

In the following sections we study spin liquid phases of 
the frustrated quantum spin model ([S]) using Schwinger 
boson mean field theory and PSG analysis. As demon- 
strated in Sec. IVIl those spin liquid phases correspond 
to FQH states of the system. 



III. SCHWINGER BOSON MEAN FIELD 
THEORY 

Schwinger boson mean field theory (SBMFT) was first 
proposed by Arovas and Auerbach^ for SU(2) invari- 
ant spin models and was further applied to frustrated 
antifcrromagnets^^ and generalized to anisotropic spin 
models.**^ Within the SBMFT spin-1/2 operators are 
represented in terms of the Schwinger boson operators 



ai-f and as follows: 

In addition, the Schwinger bosons satisfy the local con- 
straint 

4t^it + a^a^ = 1 (7) 

at every lattice site. 

To derive a mean field theory that can describe both 
ordered and disordered phases, i.e. a mean field Hamil- 
tonian that does not have any preferred direction in spin 
space, the first step is to rewrite a spin Hamiltonian 
in terms of bond operators which are bilinear forms in 
Schwinger boson operators from nearest-neighbor sites. 
Introduction of one complex Hubbard-Stratonovich field 
at each bond then provides a natural mean field decou- 
pling scheme. There are usually many possible bond op- 
erators that one can define. However, the only relevant 
ones are the bond operators that are invariant under the 
transformations of the symmetry group of the Hamilto- 
nian. The spin Hamiltonian ([5]) is invariant under spin 
rotations around the z-direction and has U{1) symme- 
try. In that case one can define four different bond op- 
erators that are invariant under symmetry operations of 
the Hamiltonian: 



F 






Aij 






Xij 






Zij 







It can be easily verified that defined bond operators com- 
mute with the operator for the total magnetization in the 
z-direction, Af = Tli^i^ ^"^^ tlms, U{\) symmetric. 
These bond operators are analogous to the bond oper- 
ators defined for the quantum XXZ model*'^ with ad- 
ditional phase factors which reflect the presence of an 
effective vector potential. 

Nonzero expectation values for these bond operators sig- 
nal the presence of various short-range correlations, that 
is, indicate the presence of short-range order. Fij and Aij 
are analogous to SU{2) ferromagnetic and antiferromag- 
netic bond operators, respectively. The bond operators 
Xij and Zij represent XY and easy axis correlations, 
respectively. Depending on the expected dominant cor- 
relations in the ground state of a given Hamiltonian, the 
Hamiltonian is further rewritten in terms of the products 
of appropriate bond operators. 

For the Hamiltonian ^ at h < he {a), where he {a) is 
a critical value of the effective Zeeman field above which 
the system is in the Mott insulating phase, the classical 
ground state has a finite total magnetization in the XY 
planed and the operators appropriate to describe the 
correlations in the ground state are Fij and Xij. The 
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Hamiltonian ([5]) can then be rewritten in the following 
form: 



(9) 



where = ni-t-+nij, = 25" is the number of bosons at each 
lattice site and is constant, and z = 4 is the number of 
nearest neighbors. This representation can then be gen- 
eralized to N Schwinger boson flavors^^'''^ and the func- 
tional integral representation of the partition function 
can be systematically expanded in powers of The 
mean field Hamiltonian corresponds to a saddle point so- 
lution of the large N action after Hubbard Stratonovich 
transformation. Since here we are interested only in the 
mean field solution, we do not use the large N language. 
We assume that the bond operators have nonzero expec- 
tation values: 



{Xij) = Xi. 



(10) 



and perform a Hartree-Fock decomposition of the Hamil- 
tonian dni): 



(11) 



where Oij G {Fij , Xij}. The mean field Hamiltonian is 
then: 



4 ^ 



(ij) i 



- 2 E (^"^t ~ ^ii) + E ^"*t + '^4 ~ ■ (12) 

i i 

Within the mean field theory the constraint on the 
Schwinger boson occupation number at site i is imposed 
on average: 



Ml 



(13) 



by introducing a chemical potential and the average 
boson number k can also be taken as a parameter4^ 
Mean field theory parameters fij and Xij are in general 
complex numbers called the mean field ansatz and can be 
found by solving self-consistent equations {Fij)MF — fij 
and {Xtj)MF = Xtj. 

Within this method, magnetic ordering in the XY 
plane is identified as Bose condensation of the Schwinger 
bosons. In the presence of the magnetic ordering the en- 
ergy excitation spectrum of Schwinger bosons is gapless, 
while a gapped excitation spectrum corresponds to in- 
compressible spin liquid phases. Further we will focus 



on the possible noncondensed ground states of the sys- 
tem which are expected to be states without any spon- 
taneously broken symmetry. As explained in one of the 
following sections, for such states we can use projective 
symmetry group analysis (PSG) to determine different 
mean field ansatze that reflect all the physical symme- 
tries of the system, and to classify possible spin liquid 
states that can appear in the phase diagram of the sys- 
tem. 



IV. MAGNETIC SYMMETRY GROUP 

Before proceeding to the PSG analysis we deflne a 
group of operators that correspond to the symmetries 
of the system and which commute with the Hamiltonian 
^ . The presence of an effective vector potential reduces 
the physical symmetries of the lattice and the Hamil- 
tonian ^ does not commute with the standard lattice 
spatial symmetry operations. However, it is still possible 
to define operators that correspond to physical symme- 
tries of the system. As will be described below, such 
operators obey the same multiplication relations as the 
original lattice spatial symmetry operators apart from 
additional phase factors and commute with the Hamilto- 
nian ([SJ. Generalized translation, rotation and reflection 
transformation operators in the presence of a uniform 
magnetic field can be specified and from such operators 
a full group can be constructed which is known as the 
magnetic symmetry group42i^"i^ 

We define elementary translation {Tx, Ty), rotation (i?) 
and reflection (r) operators, illustrated in Fig. [l] in terms 
of their commutation relations with the spin operators 
. Since we have chosen the Landau gauge, the Hamil- 
tonian ([5]) commutes with the translation operator in y 
direction and the operator Ty can be defined by 



^y^ix.y) 



C+ rp 



(14) 



It can be easily verified that [Heff,Ty]_ = 0. However, 
the operator of a pure translation in x direction does not 
commute with the Hamiltonian ^ due to the presence 
of X dependent phase factors, and we therefore define a 
generalized unitary translation operator in x direction as 
follows: 



(15) 



where the phase factors in the previous equation are cho- 
sen to ensure [ffg//, = 0. For a = p/q the Hamil- 
tonian has periodicity q in x direction and T^S^^ = 
^^x+i vY^x- Similarly one can define a unitary rotation 
operator giving rotation by 7r/2 counterclockwise around 
the origin: 



RS+ , = ^■2^o.xyg+ p 
(a:,!/) {-y,x) ' 



(16) 



where the phase factor is again chosen to obtain 
[Heff,R\_ — 0. Also, an antiunitary operator f that 
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FIG. 1; (Color online) Symmetry operations (see text). 



corresponds to reflection in the line y — x can be de- 
fined, 



i27rcKxy q 



(17) 



which in combination with the antiunitary complex con- 
jugation operator, C, forms a unitary operator and is a 
symmetry of the Hamiltonian ([5]) : 



tXS+ , = A*e-'2'ra:i=y C+ 



(18) 



where r ~ Cf , A is a complex number, and [H^f f ^t\_ — 
0. The reflection operator itself is not a symmetry of the 
Hamiltonian since it reverses chirality while the magnetic 
field explicitly breaks chirality and time reversal symme- 
try. However, after combination with complex conjuga- 
tion the appropriate sense of circulation is restored and 
operator r denotes a symmetry of the Hamiltonian. 

The magnetic symmetry group operators defined above 
have the same multiplication relations as ordinary spatial 
group operators apart from additional phase factors: 



TxR = 

TTy = 



T T 



(19) 



RTy ^ , RTx — TyR , i?/* — 1, 

Tj;T , TTj; = TyT, 



RtRt = 1 , TT = 1. 

We also note that using the Schwinger boson represen- 
tation of the spin operators (H]), the magnetic symmetry 
group operators can also be defined in terms of their com- 
mutation relations with Schwinger bosons as follows: 



TxO'(^x,y) 



(20) 



—i27ra^y rp 

— ^ '^(x+l,y)(y-'-xj 

i27ZOiaXy p 

^^{x,y)(T ^ ^[ — y,x)a-^j 



where A is a complex number, a G {t, i}, = a/2, a-f = 
—a/2 and a = a^ — a^, which is a definition convenient 
for PSG analysis. The equations (ITUl) reflect the structure 



of the magnetic space group and will play a crucial role 
in determining possible mean fleld ansatze using PSG 
analysis as described in the following section. 



V. PROJECTIVE SYMMETRY GROUP 
ANALYSIS 

In this section we discuss possible spin liquid phases 
of the Hamiltonian (|12p . Spin liquid states are spin dis- 
ordered states that typically appear in frustrated spin 
models4ii^ To study and classify possible spin liquid 
phases we use a symmetry based analysis, PSG analy- 
sis, originally introduced by Wen and collaborators 4^— 
Within this theory we show that FQH states induced by 
a lattice potential, ^'^ and with no counterpart in the con- 
tinuum, correspond to "tt flux" spin liquid states of the 
effective frustrated spin model. To the best of our knowl- 
edge the connection between the spin liquid phases of the 
effective spin Hamiltonian and FQH states for bosonic 
atoms in rotating optical lattices has not been previously 
investigated. 

As we are interested in the phases of the system that do 
not break any symmetries, we require that the Schwinger 
boson mean field theory reflects all the underlying micro- 
scopic symmetries of the spin model This will lead to 
symmetric spin liquid states. The symmetry transforma- 
tions include magnetic symmetry group transformations 
and C/(l) spin rotation symmetry. Using PSG analysis 
one can then construct all symmetry allowed mean field 
ansatze. It is important, however, to note that the micro- 
scopic symmetries of the spin model can be preserved in 
the mean field state in some indirect way. As first noted 
by Wen and collaborators^"— in the context of fermionic 
mean field theory of spin liquids, the mean field theory 
ansatz should be invariant under transformations that 
are combined physical symmetry and gauge group trans- 
formations. This idea was then generalized to Schwinger 
boson mean field theory for spin liquid states by Wang 
and Vishwanath4^i^ The transformations that leave the 
mean field ansatz invariant then form a group called pro- 
jective symmetry group (PSG). To find PSG transfor- 
mations we first note that (spin independent) local U{1) 
gauge transformations of Schwinger bosons 

a,^ ^ e"^-a,^, (21) 

leave the Hamiltonian (0) and all physical observables 
unchanged. However, the mean field ansatz is not invari- 
ant with respect to transformations (PT|) and transforms 
as follows: 



X 



J{4>i+'t>j 



i)F- 



(22) 



Since the physical spin state should be gauge invariant 
we see that that two different mean field ansatze that are 
related by a gauge transformation should, after imple- 
menting constraint p3p (that is, after projection to the 
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physical spin space) lead to the same spin state. In other 
words, different mean field ansatze may correspond to 
the same physical state and the underlying symmetries 
may not be explicitly present in the mean field ansatz. 
Presence of the local U{1) gauge transformations then 
explains why invariance of the mean field ansatz with 
respect to all microscopic symmetries is not a necessary 
condition to obtain symmetric spin liquid states. If the 
transformed (for example lattice translated) mean field 
ansatz is equivalent to U{1) gauge transformed form of 
the ansatz the same spin state will be obtained after pro- 
jection. Accordingly, invariance with respect to combined 
physical symmetry and gauge transformations is a suffi- 
cient condition to obtain symmetric states. 

In addition to PSG transformations that are related to 
physical symmetries, one finds that there are also some 
elements of PSG that are pure local U{1) gauge trans- 
formations. These transformations are not the result 
of a physical symmetry and can be associated with the 
emergent gauge group that describes obtained spin liq- 
uid phase. Such transformations form a subgroup of PSG 
which is called invariant gauge group (IGG). The micro- 
scopic theory of the spin liquid states of the system in- 
volves Schwinger boson spinous which are coupled to the 
emergent gauge field,'*°i^'^'^'* When the spinon spectrum 
is gapped and the emergent gauge field takes discrete 
values^i^ {Z2 spin liquid), or emergent gauge theory is 
U{1) gauge theory with a nontrivial topological termSS, 
the mean field solution is stable and is a good starting 
point for identifying possible spin liquid states. 

To examine the stability of a mean field solution one 
has to consider the fluctuations in Fij, Xij and Aj that 
descibe the collective excitations above the mean field 
state. Fij and Xij have the amplitude and phase fluc- 
tuations. The amplitude fluctuations have a finite en- 
ergy gap and are not important in discussing the stabil- 
ity and low energy properties of the possible spin liquid 
phases. The phase fiuctuations of the mean field ansatz 
are gapped in the presence of a Chern-Simons term or 
spinon- pair condensation ( Anderson- Higgs mechanism). 
Here the latter is achieved if both Fij and Xij are nonzero 
and in that case the emergent gauge theory is Z2 (Ising) 
gauge theory. If Fij or Xij is zero the emergent gauge 
theory is t7(l) gauge theory. However, a Chern-Simons 
term present in the emergent gauge theory due to an 
effective uniform magnetic field gives the U{1) gauge bo- 
son a nonzero gap. Gapped gauge bosons can mediate 
only short range interactions between spinous and stable 
spin liquid states always contain spinons with only short 
ranged interactions between spinons. 
Further we will assume Z2 spin liquid states, that is, we 
will initially assume that both Fij and Xij are nonzero. 
It is then clear that the only two elements of the IGG are 
the identity operation 1 and the IGG generator —1. In 
other words, the only two transformations that leave the 
mean field ansatz invariant and are pure local U{1) trans- 
formations (I^T|) are bi — >■ e'^'&i with (pi — or (fii = tt. The 
spin U{1) rotation symmetry is already realized by con- 



sidering mean field ansatz, Xij and F^j, of the form ([5]) 
which is explicitly invariant under spin rotation around 
z direction. The operations that we further need to con- 
sider are magnetic symmetry group operations. 

As discussed before we will consider symmetric spin liq- 
uid states that are invariant under PSG transformations. 
A PSG transformation will be a combined transforma- 
tion of magnetic symmetry group and gauge transforma- 
tions. However, it is not possible to combine a magnetic 
symmetry transformation with any gauge group trans- 
formation since algebraic relations (|19p between mag- 
netic symmetry group elements constrain the possible 
choices of gauge transformations. For each transfor- 
mation T {T € {Tx,Ty, R,t}) in the magnetic symme- 
try group there is a gauge group transformation Gt 
(Gt e {GT^,GTy,GR,Gr}) such that the mean field 
ansatz is invariant under combined transformation GtT . 
The gauge transformations can be represented as follows: 

Gt ■■ bpa ^ e'^^^^hp,. (23) 

The transformations Gt can be found by consid- 
ering algebraic constraints (|19p (with a ^ 
when the magnetic symmetry group transformations 
are applied to the Schwinger boson operators bi^). 
For example, since e'^''°'''T-^TyT^T-^ = 1, then 

e^2.<.^ (Gt^T.)-^ GT.TyGT^T, {GT.Tyy^ has to be 
equivalent to an identity operator, that is, it has to be 
an element of IGG (1 or —1). Using relation 

Y-'GxYbp, {Y-'GxYy' = e'^-I^^lfo,-,, (24) 

where here X, Y E {Tx,Ty}, we obtain the following 
equation 

[T, (r}H-0T„ [Tx m+^T. [T^Ty-^ (r)]-0T„ M^PiTT, 

(25) 

where pi £ {0, 1} corresponds to IGG elements 1 and 
-1 {e'P''^ with pi = or 1). 

Similar equations can be obtained by considering other 
algebraic relations in the equations (fT9)) . A detailed 
derivation of those equations is described in Appendix |X] 
and the most general solution that satisfies the obtained 
equations is: 

<^T. [x,y] = 0, (26) 

(j>r [x, y] = (l)r [0, 0] -I- PlTTXy, 
(l}R[x,y] = (l)R[0,Q]+piTrxy, 

where pi £ {0, 1}. If the nearest neighbour amplitudes 
Xij are nonzero then (/lij [0, 0] = 0r [0,0] = p27r/2 with 
P2 G {0,1}. Possible mean field ansatze can then be 
obtained by requiring invariance of the mean field ansatz 
with respect to PSG transformations: 

{GtTF,, (GtT)-') = (F,,) = (27) 
{GtTX,, (GtT)-') = (X„-) = X,,, 
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where i and j are nearest neighbours, T € {T^^ Ty, R, r} 
and fij and Xij are both real and nonzero in general 
(Appendix |A|. From the conditions ([27| we obtain four 
possible mean field ansatze: 



C : f{xi,yi) 

^ ■ f{xi,yi) 
~^(xi,yi) 

where = /*,, 



f(xi,yi){xi,yi + l) - fi (28) 

.f{xi,yi)(xi,yi + l) — /j 
*^(xi,yi)(a:^,yj + l) 

( — 1)^" X(.j;i,yi){xi + l,yi) — X, 
(^1)^* /(a;.,i/,)(2:i + l/yi) = /, 



a; 



{xi + l,yi 



G M and above equa- 



tions are valid for all lattice sites i = {xi,yi). The mean 
field ansatze above, illustrated in Fig. [5J then give local 
minima of the mean field Hamiltonian. The ansatze A 
and B can be distinguished from C and D by the emer- 
gent gauge invariant flux through an elementary square 
plaquette defined as the gauge invariant phase of the fol- 
lowing products^ : 



fijfjkfklfu = I/I , 
•^ijXjk-^kl^li — l-^l 



(29) 



with i, j, k, I being sites of a square lattice plaquette. For 
the ansatze A and B the phase $ = (zero-flux states) 
and for C and D the phase $ = tt (tt-Aux states) for 
all square lattice plaquettes. Therefore these two sets of 
states are clearly not gauge equivalent meanfield states. 
Also, A and B differ by the symmetry of the Schwinger 
boson pairing order parameter Xij. Unlike for ansatz A, 
for ansatz B the order parameter Xij changes sign un- 
der a 7r/2 rotation (d wave like state). For 7r-flux states, 
C and D, the unit cell for Schwinger bosons is doubled 
due to the emergent tt flux through each unit square pla- 
quette of the lattice and the translation symmetry is not 
explicit in the mean field Hamiltonian. Similarly as for 
the ansatze A and B, ansatze C and D differ by the 
symmetry of the order parameter Xij. In the following 
section we further study possible noncondensed phases 
arising from the ansatze described previously and relate 
these states to bosonic fractional quantum Hall states for 
bosonic atoms in rotating optical latticesJ^ 



VI. BOSONIC FRACTIONAL QUANTUM HALL 
STATES 



To further study noncondensed phases of the mean 
field Hamiltonian we take into account the local con- 
straint ([7]) on Schwinger boson number at each site 



B 


i 1 


t < 





FIG. 2: The zero-flux (A and B) and the 7r-flux (C and D) 
mean field ansatze. The solid black, solid gray and dashed 
lines indicate bonds with fij = / and Xij = x, fij — f and 
Xij = —X and fij = — / and Xij — ~x, respectively. Here 
fij — f'ji, Xij — Xji and both / and x are real. 



exactly;^"— unlike in the conventional Schwinger boson 
mean field theory where the constraint is imposed on av- 
erage as in the equation p^ . The local constraint ([7]) 
introduces strong correlations between spin-up and spin- 
down bosons. First we note that the constraint implies 
that the number of up/down bosons at each site can be 
only zero or one. Thus we require that the Schwinger 
boson operators satisfy the following anticommutation 
relations on the same site: 



{c 



'}+ = > {a»,T,a[^'| = S„ 



where cr, cr' € {tji}? a-nd commutation relations 



a 



= 0, 



(30) 



(31) 



for the lattice sites i and j when i =/= j. The requirement 
above corresponds to a hard-core limit for each bosonic 
species and excludes states with more than one boson of 
the same species at the same lattice site. However, to 
take the constraint ([7]) into account exactly we also need 
to exclude the states with one spin-up and one spin-down 
Schwinger boson at the same lattice site. This is achieved 
by the Gutzwiller projection: 



(32) 



where n^^ and n^^ are the number operators at a site i 
for spin-up and spin-down Schwinger bosons respectively, 
which for hard-core bosons have eigenvalues or 1. The 
ground state wave function is then 



ip = Palp 



MF, 



(33) 
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where ipMF is the mean field Hamiltonian ground state itly it is convenient to rewrite the mean field Hamiltonian 
wave function. To apply the Gutzwiller projection explic- p2p in the following form: 



"I E^^^- {^'^''^' + h^tPl) + e-'^-/' {h,,pl+Puhl,) } , 

(ij) 

I 



where 

ho = j^{\fij\^ + + + 

(ij) (ij) i 

(34) 

with Ns being the number of lattice sites, and where we 
have introduced the notation pi^^ = a^j, and hi^ = a\^. 
We have also taken into account that fij and Xij are real 
for all mean field ansatze (pS)) . The operator creates a 
spin-up bosonic hole and the operator creates a spin- 
down bosonic particle at the site i. The constraint ([7]) 
can then be rewritten in the form: 

p\iPii = ^It^^t. (35) 

where we have used the anticommutation relations pop . 
The local constraint is equivalent to saying that 

wherever there is a spin-down boson there is also a spin- 
up hole. The Gutzwiller projected wave function with 
N]^ spin-down bosons and TV — iV^ spin-up bosons can 
then be obtained in the following way. Starting from 
the state with N spin-up bosons, first spin-up bosons 
are annihilated and then spin-down bosons created 
at the coordinates of the spin-up holes. Therefore the 
wave function in the coordinate representation will be 
a function of the coordinates of spin-down bosons (or 
equivalently of spin-up bosonic holes) only. 

Before examining further the effect of the Gutzwiller 
projection we note that the Hamiltonian ([M)) describes 
correlated hopping of particle-hole pairs. We are thus led 
to introduce the following hard core boson operators: 



bi± ^ —7= {Pii ± /i»t) ■ (36) 



It can be easily verified that these operators satisfy mixed 
commutation - anticommutation relations of the form 
(150)) and (PT|) . In terms of operators ([55]) the mean field 



Hamiltonian can be written as: 

Hmf = + ^ E + ^1-^^-) (37) 

i 

where tfj ~ fij ± Xij. Also, the constraint ([35]) implies 
the following condition: 

blh+ = blh^. (38) 

This can be seen from the equation fej^6i+ — b\_bi^ = 

Pi-l^it + ^l^Pii- When the constraint (|55)) is taken into 
account explicitly then the only allowed states are the 
states with zero or one particle-hole pairs at a lattice site. 
Within this restricted single site basis terms p^ft-j^ ^^nd 

hl^Pil can be set to zero and condition psp is obtained. 
The ground state wave function can then be obtained 
through the following projection: 

N N 

^{{n}) = {0\Y[h,^p,i\MF) = {0\l[b,+b,^\MF}, 

i=l i=l 

(39) 

where |0) represents the vacuum of spin-down bosons 
and spin-up holes, (or equivalently vacuum of 64- and 
6_ bosons), Ti are coordinates of N bosonic atoms in the 
lattice, and \MF) is a mean field state of bosons 6+ and 

The effective fiux per plaquette of the lattice for 6+ /5_ 
bosons is ot^fj = Oieff/'^, with Uefj being the effective 
flux for the atoms. For the mean field ansatze A and 
B the effective flux afff = ce/2, and for the mean held 
anzatze C and D ct^ff = Q;/2 -I- 1/2. Since the number of 
atoms is the same as the number of 6+ (6_) bosons and 
oi^ff = O'-eff/'^ the wavefunction can be written as: 

^A{n})^W.,=2A{m\ (40) 

where ■0^, ({ri}) is the wavefunction for b+/b^ bosons 
at filling fractions v' = n+/a^fr = n^/a^rr = 2v = 
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2n/aeff: with n_|_ /n_ being the average density of 6+ /6_ 
bosons per lattice site. Thus at an atom filhng frac- 
tion V, the average density of atoms per lattice site is 
n = av' /2 for the mean field ansatze A and B, and 
n = (1/2 ± a/2) v' for the mean field ansatze C and _D, 
where we have taken into account the symmetries under 
cf^j j — ^ 1 — cf^j j and n — > 1 — ri (particle-hole symmetry 
on the lattice). 

We further argue that the FQH states of bosonic atoms 
in a rotating optical lattice correspond to the incompress- 
ible states for 6+ — &_ pairs. Without the constraint, two 
species of hard-core bosons, &+ and 6_, are independent 
of each other and the ground state wavefunction of the 
mean-field Hamiltonian (|37p can be written as a product 
of the ground state wavefunctions for 6+ and &_ bosons. 
For example, at an effective filling fraction v'^^^ = 1 the 
ground state for each bosonic species (6+ and &_) is pre- 
dicted to be an incompressible state that, in the con- 
tinuum limit, can be described by the bosonic Pfaffian 
stateji The projected wavefunction (|40|) then corresponds 
to the z^e// = 1/2 FQH states for 6+ — &_ pairs (or equiv- 
alently for particle-hole pairs). 

At v'pjj — 1, that is v^ff — 1/2, possible incompressible 
ground states are at n = a/2 (mean field ansatze A and 
B) and n = l/2±a/2 (mean field ansatze C and D). 
The states with n = l/2±a/2 are lattice induced FQH 
states, with no counterpart in the continuum, and within 
presented mean field theory, correspond to "tt flux" spin 
liquid states of the effective spin model (mean field an- 
satze C and D) for which the hard-core bosons 6+ and 
6_ see an effective flux oc^^^ = 1/2 + a/2. Such incom- 
pressible quantum Hall states have previously been found 
using the composite fermion approach and exact diago- 
nalization studies Several other lattice induced incom- 
pressible states predicted within the composite fermion 
approach can also be related to incompressible states for 
6+ — 6_ pairs, for example states at i^'^^ff = 2/3. 
Finally, we note that the projected wavefunctions are also 
classified by PSG and that our results provide a guide- 
line for a further search of possible ground states in the 
projected wavefunction space. We thus propose to study 
the effective spin model by Gutzwiller projected wave- 
function variational approach^"— as a possible direction 
of future research. 



VII. CONCLUSIONS 

We have studied possible incompressible ground states 
for bosonic atoms in a rotating square optical lattice. In 
the hard-core boson limit the Bose-Hubbard Hamiltonian 
for the system can be mapped to a frustrated easy-plane 
magnet. Using Schwinger boson mean field theory and 
projective symmetry group analysis we have classified 
possible stable gapped spin liquid states of the effective 
spin model and related those states to FQH states for 
bosonic atoms in rotating optical lattices. In particular, 
we have found that, within the mean field theory devel- 



oped here, the FQH states induced by lattice potential 
and with no counterpart in the continuum limit, corre- 
spond to "tt flux" spin liquid states of the effective spin 
model. Since there are also competing condensed states 
on the latticei^ii^'^ it is important to test the predic- 
tions of the mean field theory presented here. Finally, 
further work is necessary for more direct comparison of 
the obtained results with known results for FQH states 
in rotating optical lattices. We hope the Gutzwiller pro- 
jected wavefunction variational approach will be further 
developed in future research. 
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Appendix A: Solution to PSG equations 

In this appendix we derive equations for allowed gauge 
transformations Gt, defined by equation ([23|) in section 

m 

Gt ■■ bp, ^ e^t'^'^'^hp,. (Al) 

As explained in section |V] the mean field ansatz should 
be invariant under combined transformations GtT where 
T are magnetic symmetry group transformations. The 
transformations Gt can be found by considering alge- 
braic constraints (|19p following the procedure explained 
in the example in the main text. We first describe in 
more details derivation of the equation ( [25|) and then ap- 
ply the same procedure to all other algebraic constraints 
given by equations (|19p . 
From 

ei'^^°''T-%T^T-^ = 1 (A2) 

it follows 

e^^-"'{GT^T^r^GTTyGT^T^{GTTy)-^=e'P'^ (A3) 

where pi E {0, 1} corresponds to IGG elements 1 and —1. 
The relation (jA3l) can be rewritten as 

X '-'T^^^^x '^Ty-lxJ^x ^V-Lx-Ly lylx '-'T^J^xJ^y <^Ty 
^^^p,n^-^2na^^ (A4) 

Since T^T"! = 1, and T'^TyT^Ty^ = g-^^^"- and 

Y~^G,YW, (r-iG,r)"' = e'-^-I^^lfo,-,, (A5) 
we obtain equation ([25]): 

-<j,T^ [T, (rO]+0T„ [Tx (r)]+</.T. [T.T-i (f)]-(/.T„ [r1 =pi^, 

(A6) 
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with pi G {0, 1}. We also note that, since (j)T are phases, 
all equations for phases in this section will be true modulo 
27r. If a gauge transformation bpcr — > e^'^'^'^^^bp^j is applied 
then phases in (|A1[) change as follows: 



0T {r) ^ <t>T (r) + c^G {r) ~ <I>G [T~' (r)] 



(A7) 



Using this gauge freedom the mean field ansatz can be 
made explicitly invariant under translation in x direc- 
tion, that is we can assume 

{r) = 0. (A8) 

Then the equation (|A6P becomes: 

</>T„ [2: + 1,2/]- <l>v [x, y] = PiTT. (A9) 

The solution of the equation above is then: 

(t^Ty [x, y] = 0T„ [0, y] + piTTx. (AlO) 

We further assume that using gauge freedom we can set 
(j)Ty [0,y] — while keeping (r) — 0. These two equa- 
tions can be simultaneously satisfied for open boundary 
condition. We further assume open boundary condition 
in which case equation (jAlOp becomes: 



(f>Ty [X,y] =PlTTX. 



(All) 



Following the same procedure, from the algebraic con- 
straints T-^RT-^R-^ = 1 and T'^RT^R-^ = 1 we ob- 

X y y ^ 

tain: 



[x + l,y]- (f>R[x,y]=p'2T^ + PiT^y, (A12) 



and 



ipR [x, y + l\-(j)R [x, y] = pgTT + pi-Kx, (A13) 

with p'21 P3 G {0, 1}. The solution of the above equations 
is: 

0i? \x, y\ = (\>R [0, 0] -f p'^nx + pgTry -|- piTrxy, (A14) 

Equivalently from rTj, = T^^t and tT^ = TyT (or fTy = 
TxT and e'^'^'^'^'^fTx — Tyf) we obtain foUowing equa- 
tions: 



and 



<j>r[x + l,y]-(j)r[x,y]=p'^TT+piTry (A15) 



(j>r [x, y+l]~(j>r [x, y] = PgTT + piTTX, (Al6) 



with P4, P5 e {0; 1}. The solution of the equations (jAlSP 
and (|A16|1 is: 

0r [X, y] = cj)r [0, 0] -I- p'^TTX + p'^Tiy -f PlTTXy. (A17) 



If we further apply the gauge transformation <j)Q [x,y\ = 
TTx the phases (j)T will change as follows: 

<I>T^ [X, y] <I>T^ [X, y] + TT = TT, (A18) 

0T„[a;,y] (j>Ty [x,y] ^ piirx, 

<I>R [x, y] <I>R [x, y] + IT {x - y) = (f>]i [0, 0] -I- {p'2 + 1) nx 
+ {p'3-^)^y + Pi'^xy, 

(pr [X, y] (j>r [x, y] + TT {x - y) ^ cl)R [0, 0] + (P4 + 1) TTX 

+ {p'5 — 1) Try + pinxy. 

This gauge transformation does not change (j)T {r), nor 
does it really change (pT^ (r) since a site-independent con- 
stant TT can be added to any (f>T (r) because of the IGG 
structure. However, when the transformation is applied 
P2,P4 ^P2 + l>-P4 + l: andp3,P5 ^^3-1,^5-1. Since 
Pi = and p'i — I are gauge equivalent V i G {2, 3, 4, 5} 
we can always assume p'2 ~ P-} = P4 = p'^ ~ 0- 
Then we can write: 



'/'T, [x, y] 
4>Ty [x, y] 
pR [x, y] 
0T [x, y] 



0, 

<pR [0,0] -f piTTxy, 

(j)r [0, 0] -I- PlTTXy. 



(A19) 



The allowed mean field ansatze can then be obtained by 
requireing invariance of the mean field ansatz with re- 
spect to PSG transformations: 



(A20) 



where i and j are nearest neighbours and 
T e{T,,Ty,R,T}. 

Considering further following equations for 
the (0,0) (-1,0) bond: /(o,o)(-i,o) 

and 



0)(1,0) 



/(-1,0)(0,0) = /(0,0)(1,0) = /(0,0)(-l,0) '^'^ ^"^"^ ^^^^ 

/(o,o)(-i,o) = /(o,o)(-i,o) = /■ III otli^^i' words /(o,o)(-i,o) 
has to be real. Also, a global phase rotation, that is, a 
gauge transformation (pa (r) = (f) does not change any 
PSG elements (equation (jA7| ). This can be used to 



fix one of the Xij presented to be real and positive, for 

example, 2:(o^o)(o,i) = ^(o,o)(o,i) = ^■ 

If we assume that the nearest neighbour amplitudes 

Xij are nonzero then [0, 0] can be found from the 



following equation for the (0, 0) 

a^(o,o)(~i,o) 
which in 



(-1,0) bond: 



with the equations 
2:(o,o)(i,o) (due to the 



combination 

X{-ifi){ofi) 



X(0,0)(-lfi) = 

symmetry Xij — Xji and translational invariance 
of direction) gives ^(pR [0, 0] = 0, that 

is, 0fl[O,O] = P27r/2 with p2 € {0,1}. Equiva- 



lently from X 



(0,0)(0,1) 



-j((/.h[O,O]-I-0r[1,O]) 



2;(0,0)(1,0) 



e-*('^-[°'°I+'^-[-i'°l)a;(o,o)(-i.o) we obtain 0r[O,O] = 
(pn [0,0]. The algebraic constraints i?"* = 1, RtRt = 1 
and rr = 1 do not further restrict possible values 
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of Pi and p2- The mean field ansatze ^7} can then (pu [0,0] = (j)r [0,0] =P2Tt/2. 
be obtained from equations (jA19p and (|A20|) with 
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